Flux Qubits and Readout Device with Two Independent Flux Lines 
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We report measurements on two superconducting flux qubits coupled to a readout Supercon- 
ducting QUantum Interference Device (SQUID) . Two on-chip flux bias lines allow independent flux 
control of any two of the three elements, as illustrated by a two-dimensional qubit flux map. The 
application of microwaves yields a frequency-flux dispersion curve for 1- and 2-photon driving of 
the single-qubit excited state, and coherent manipulation of the single-qubit state results in Rabi 
oscillations and Ramsey fringes. This architecture should be scalable to many qubits and SQUIDs 
on a single chip. 

PACS numbers: 03.67.Lx, 85.25.Cp, 85.25.Dq 



Superconducting quantum bits (qubits) based on 
charge 0, Q , magnetic flux [3, |3 , and phase difference 
across a Josephson junction [j, |a| are attractive candi- 
dates for the basis of a quantum computer because of 
their inherent scalabihty using established thin-film fab- 
rication techniques. Advantages of the flux qubit include 
its immunity to the ubiquitous charge noise in the sub- 
strate and that it can be conflgured with no direct elec- 
trical connections. One type of flux qubit consists of 
a superconducting loop interrupted by three Josephson 
junctions with critical currents Iq, Iq, and alo (a < 1) 7]. 
When the applied flux bias $q is at a degeneracy point 
(n-|-l/2)$o (n is an integer such that |$q— n$o| < "^o/S; 
$0 = h/2e is the flux quantum), a screening supercur- 
rent Jq can flow in either direction around the loop. 
The ground and first excited states of the qubit corre- 
spond to symmetric and antisymmetric superpositions 
of the two current states and are separated by an en- 
ergy A. Here, A//i is the tunnel frequency between the 
current states, typically a few GHz. When is away 
from a degeneracy point, the energy difference between 
the two superposed states is = (A^ -I- e^)^/^, where 
e = 2Jq[$q - {n+ l/2)<i>o]. The state of the qubit is 
measured by coupling the screening flux generated by Jq 
to a hysteretic dc superconducting quantum interference 
device (SQUID). This flux determines the bias current at 
which the SQUID switches out of the zero-voltage state. 

In addition to the development of scalable interqubit 
couplings a prerequisite for scaling to a system of 
many qubits is that the attendant readout, filtering, and 
bias circuitry also scale. A particular challenge is that 
the flux bias must be settable for each element individ- 
ually. This mandates the use of on-chip flux-bias lines 
in an arrangement that enables one to apply a combi- 
nation of currents to address any given qubit or SQUID 
while maintaining all other flux biases at constant val- 
ues. Furthermore, the bias currents required to change 
the flux over (say) ±l$o should not be so large that it be- 
comes impractical to deliver them to a chip cooled to mil- 
likelvin temperatures. This requirement establishes min- 
imum self-inductances of the qubit and readout SQUID 



that are substantially larger than values used previously 
in 3-junction qubits, which have relied on external coils 
to generate large magnetic fields 0, 0, ^3 ■ the same 
time, the mutual inductance between the on-chip flux 
lines and the qubit must be sufficiently small for the noise 
generated by the circuitry supplying the flux bias current 
not to be the limiting source of decoherence ■ 

In this Letter, we report measurements on two qubits 
and a readout SQUID that meet these criteria. We illus- 
trate the orthogonalization of the applied fluxes by means 
of a two-dimensional flux map. We report spectroscopy 
on one of the two qubits that matches the expected dis- 
persion for a flux qubit. In addition, we observe spurious 
resonances, some of which may be related to defects in 
the tunnel barriers of the qubit. We perform coherent 
manipulation of the single-qubit state resulting in Rabi 
oscillations and Ramsey fringes. In prior measurements 
of coherent oscillations 0,^3' readout SQUID was 
connected directly to the flux qubit and thus detected a 
combination of flux and phase changes. In contrast, our 
SQUID is electrically isolated from the qubit and detects 
only flux changes. 

Figure 1 shows the device layout. The readout SQUID 
has a calculated inductance Ls — 358 pH, and each of 
the two qubits it encloses has a calculated inductance 
Lq = 143 pH. The calculated mutual inductance between 
each qubit and the SQUID is 61 pH. One pair of series- 
connected flux bias lines is arranged near the top of the 
SQUID and a second near the bottom; thus, flux in qubit 
1 (2) is supplied predominantly by the lower (upper) flux 
lines. The mutual inductance between each qubit and its 
associated flux lines was designed to be 4-5 pH, enabling 
us to apply ~ 2$o with a current of 1 mA. This criterion 
dictated the relatively large qubit self-inductances com- 
pared with those in previous experiments 0, 0, 0, 0| . 

We fabricated the device on an oxidized Si substrate 
using electron-beam lithography and double-angle evap- 
oration to form the Al-AlOx-Al tunnel junctions. The 
Al lines for the qubit and SQUID loops were 1 /im wide, 
and those for the flux bias 10 /im wide. Each SQUID 
junction was 175 x 200 nm^ with a critical current of 
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FIG. 1: (a) Chip layout. Dark gray represents Al traces, 
light gray AuCu traces. Pads near upper edge of chip provide 
two independent flux lines; wirebonded Al jumpers couple 
left and right halves. Pads near lower edge of chip supply 
current pulses to the readout SQUID and sense any resulting 
voltage, (b) Photograph of center region of completed device. 
Segments of flux lines are visible to left and right of SQUID, 
which surrounds the two qubits. 



220 nA. For qubits 1 and 2, the larger junctions had ar- 
eas of 250 X 250 nm^ and 180 x 200 nm^, approximate 
critical currents Iq (scaled from SQUID junction areas) 
of 390 nA and 230 nA, and a-values (based on junction 
areas) of 0.49 and 0.68, respectively. The different junc- 
tion parameters were chosen to increase the probability 
of obtaining one viable qubit; in fact, qubit 2 displayed 
good characteristics. A 42 nm-thick AuCu film deposited 
and patterned prior to the Al deposition provided quasi- 
particle traps near the junctions T^l, 100 fl shunts on 
each of the four flux lines, and 500 fl and 1275 fl series 
resistors on each end of the pulse lines and sense lines, 
respectively. 

To eliminate external magnetic field fluctuations we en- 
closed the chip in a 6 X 16 X 22 mm'^ cavity machined into a 
copper block and plated with Pb. This superconducting 
cavity stabilizes the field, enabling us to acquire data for 
periods up to 48 hours. A 1 mm-diameter superconduct- 
ing loop ^ 3 mm above the chip supplied microwave flux. 
The sample holder was attached to the mixing chamber 
of a dilution refrigerator at 50 mK. All electrical leads 
to the experiment were heavily filtered at several differ- 
ent temperatures with a combination of lumped circuit 
and copper powder low-pass filters jl5j . Measurements 
of the state of a qubit as a function of flux, microwave 
power and frequency were made by pulsing the current 
in the SQUID and detecting whether or not it switched 
out of the zero-voltage state by means of a low-noise, 
room-temperature amplifier. The flux bias currents were 
supplied by highly stable potentiometers, controlled by a 
computer over a fiber-optic link [l^ . 

There are three applied fluxes that determine the state 



of the system: the SQUID flux $5, and the qubit fluxes 
^Qi and $Q2. Given the two flux lines, we can set any 
two fluxes arbitrarily; the third is fully constrained. This 
can be expressed succinctly by the matrix equation 
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where the Af^ are mutual inductances between the vari- 
ous fiux lines and loops, the <I>° account for static back- 
ground fields, and the li are currents in the flux lines. 

We turn now to our experimental results. Figure 2(a) 
shows the switching probability of the SQUID versus the 
amplitude of the current pulses applied to it for two dif- 
ferent values of $q2 at constant $s. Because we hold 
$s constant, we can measure the displacement of the 
two curves for constant sensitivity in the SQUID. The 
measurement fidelity, which is the difference between the 
switching probabilities, has a maximum of about 60%. 
Figure 2(b) shows /f"^" vs. $s; -f™^" is the pulse am- 
plitude for which the switching probability is 50%. The 
effects of the changing flux in the qubits are small on this 
scale. Figure 2(c) shows vs. $qi for constant $s; 
as ^Qi is varied, a flux of approximately the same mag- 
nitude and opposite sign is applied to qubit 2. Thus, as 
<i>Qi is increased at constant ^>s, abruptly increases 
when qubit 1 flips state and decreases when qubit 2 flips 
state. 

To determine the parameters in Eq. J^l, we sampled 
at - 20,000 different settings of h and h- We flt 
these data to a parametric model of the response of the 
SQUID //**($) to the total SQUID flux $t- To describe 
the SQUID modulation, we use the ad hoc expression 
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*($) = Vai cos 27ri h V 6* sin 27ri h d, (2) 

j=l 1=2 



where the Oi, 6^, and d are fit parameters. The total flux 
coupled to the SQUID is well approximated by 

*T = + A$QijQi($Qi) + A$Q2jQ2(«'Q2), (3) 

which neglects the self-screening of the SQUID. Here, the 
A$Qi = JgiMqis are the amplitudes of qubit screening 
flux changes near the degeneracy point, referred to the 
SQUID, and the qubit circulating currents are described 
by dimensionless periodic functions of unit amplitude 

JQ»(*) = {sin[27rgQ,($-n$o)/*o]}/sin7r5Qi, (4) 

with fitting parameters ggi. The jQi{^) are discontinu- 
ous when n increments. 

Fitting the /f ^"(/i, J2) data, we obtain the following 
parameters for //**($t): Mpigi = 3.96 pH, Mf2Qi = 
-0.77 pH, Mfiq2 = 1.37 pH, Mf2Q2 = -3.31 pH, 
Mfis = 8.30 pH, Mf2S = -6.30 pH, A$qi = 4.39 
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FIG. 2: (a) SQUID switching probability vs. amplitude of 
bias current pulse near qubit 2 transition. The two curves 
represent the states corresponding to $q2 = 0.48'I>o (red) 
and 3><32 ~ 0.52"I>o (blue); $s is held constant. Each curve 
contains 100 points averaged 8,000 times, (b) 7^''^° vs. <E>s. 
Each period of oscillation contains ~ 5, 000 flux values, and 
each switching current is averaged 8, 000 times, (c) Depen- 
dence of /g"'^ on "I>Qi for constant "I>s. (d) Qubit flux map. 



m$o, A$Q2 = 4.04 m$o, 9qi = 0.64, gQ2 = 0.51. This 
fit captures all the essential features of the data, having 
a root-mean-square residual of 1.8 nA, that is, less than 
0.5% of the maximum value of If 

In Fig. 2(d), we plot I^"'^" - //'*($5) versus the cur- 
rents Ii and I2; lines of constant $s, ^qi, and $q2 
are indicated. This qubit flux map displays only the 
flux contributions of the two qubits. For example, in- 
side the square, where $qi « $0 and <i>Q2 ~ — $o/2, 
j50% exhibits an abrupt step across the $q2 line, as Jq2 
changes direction discontinuously. However, there is no 
step across the $qi line, as Jqi is continuous. Further- 
more, one can scan the switching current along a line 
where $s is held constant at an arbitrary value, to ob- 
serve the effects of flux only on the qubits. Finally, spe- 
cial points where the two qubits are each at a degeneracy 
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FIG. 3: Spectroscopy of qubit 2. Enhancement and suppres- 
sion of if ° is shown as a function of "1>q2 and fm relative 
to measurements in the absence of microwaves. Dashed lines 
indicate fit to hyperbolic dispersion for 1- and 2-photon qubit 
excitations. The 2-photon fit is one-half the frequency of the 
1-photon fit. Inset containing ~ 23, 000 points is at higher 
resolution. 



point are indicated by the circle. 

The remaining discussion of the experiments is con- 
cerned with the effects of an applied microwave flux on 
qubit 2. Photons of energy hfm = v drive transitions be- 
tween the two qubit states, producing peaks and dips on 
the qubit transitions. By measuring J^"'^" as a function 
of ^Q2 and /m, we determined the dispersion relation 
shown in Fig. 3. This measurement contains a total 
of ~ 75, ODD points and took 48 hours to acquire, thus 
demonstrating the excellent flux stability in our system. 
The dispersion is well described by the hyperbolic rela- 
tion V = {(? ^ A2)i/2^ with A//i = (3.99 ± 0.05) GHz 
and {l/h)de/d^Q2 = (896 ± 5) MHz/m$o. Part of the 
spectrum corresponding to 2-photon transitions is also 
shown. 

Figure 3 shows two types of deviation from ideal be- 
havior, more evident in the inset which was measured 4 
weeks before the full spectrum. First, there are sharp 
suppressions of the critical current which are indepen- 
dent of $Q2 and hence occur at particular constant val- 
ues of /m, for example, at 7.4 GHz, that we believe arise 
from electromagnetic modes which couple to the qubit. 
Second, we see disruptions of the dispersion curve that 
are suggestive of coupling between the qubit and other 
two-state systems. One instance of this second anomaly, 
originally centered near 7 GHz in the inset and shifted 
to 6.5 GHz in the full spectrum, is remarkably similar to 
those reported for phase qubits, and may be of the same 
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FIG. 4: Coherent manipulation of qubit state, (a) Rabi 
oscillations, scaled to measured SQUID fidelity, as a function 
of width of 10.0 GHz microwave pulses, (b) Rabi frequency vs. 
10.0 GHz pulse amplitude; line is least squares fit to the data, 
(c) Ramsey fringes for qubit splitting of 9.95 GHz, microwave 
frequency of 10.095 GHz. (d) Ramsey fringe frequency vs. 
microwave frequency. Lines with slopes ±1 are fits to data. 



origin '171. 

We performed coherent manipulation of the qubit state 
by varying the duration of a resonant microwave pulse 
for fixed frequency and amplitude. Upon applying fixed- 
amplitude bias current pulses, we observed Rabi oscilla- 
tions in the switching probability of the SQUID as a func- 
tion of the microwave pulse width at approximately 100 
frequencies ranging from 4.37 to 13.86 GHz. Generally 
speaking, the amplitude of the oscillations is higher and 
their decay time is longer for qubit flux bias points away 
from the spurious splittings shown in the spectroscopy. 
An example is shown in Fig. 4(a), where each of the 
195 points was averaged 20, 000 times. We reference the 
amplitude of these oscillations to the qubit by scaling 
to the SQUID measurement fidelity. We observed Rabi 
oscillations with the longest decay time when we oper- 
ated at a bias current of 296 nA, where the fidelity was 
14.4% for this particular value of $5 and the SQUID 
switched infrequently [Fig. 2(a)]. We believe that this 
improvement is related to noise currents in the SQUID 
loop which couple to the qubit through Mq2s and are 
produced by the quasiparticles generated in the SQUID 
loop during each switching event: a lower SQUID switch- 
ing probability results in fewer quasiparticles averaged 
over time 13 ■ Switching probability data in Fig. 4(a), 
scaled to the measurement fidelity of the SQUID, fit well 
to a damped sinusoid with a visibility of 63% and a de- 
cay time T/jabj = 78 ns. Fits to similar measurements for 
different microwave pulse amplitudes show that the Rabi 
frequency scales linearly with the microwave amplitude 
[Fig. 4(b)], as expected for coherent driving [l^ . 

We measured the decay of a resonance peak and ob- 



tained the relaxation time tr = 281 ns at 10.0 GHz from 
an exponential fit. We measured the dephasing time 
from Ramsey fringes We first applied a 7r/2 pulse to 
tip the qubit state vector into the equatorial plane, where 
it dephased. To calibrate the it/ 2 pulses, we chose a mi- 
crowave pulse width based on a Rabi oscillation measure- 
ment made at the same operating point. After a variable 
time delay r, we applied a second 7r/2 pulse followed by 
a measurement of the SQUID switching probability for a 
fixed bias current pulse amplitude. When we chose a mi- 
crowave frequency off resonance we observed damped os- 
cillations of the SQUID switching probability [Fig. 4(c)] 
with a fit decay time of the oscillations — 6.6 ns. 
Figure 4(d) shows the Ramsey fringe frequency for 100 
values of the microwave frequency, together with the fit 
to two lines with slopes of magnitude unity, as expected 
for coherently driven Ramsey fringes. 

In conclusion, we have demonstrated on-chip fiux bias 
lines which allow us to vary the flux applied to two of 
the three devices independently. The mutual inductance 
between the flux lines and the qubits is weak enough that 
the characteristic impedance of the flux lines (say 50 fi) 
should not limit our coherence times Although in 

this Letter we have concentrated on the quantum coher- 
ent properties of a single flux qubit, we note that the 
two flux qubits in this design in principle could be cou- 
pled controUably using the circulating currents in the dc 
SQUID 0. This system of qubits, SQUID and flux lines 
should be readily scalable. 

We thank M.H. Devoret, D. Esteve, C.J.P.M. Har- 
mans, J.M. Martinis, R. McDermott, J.E. Mooij, R.J. 
Schoelkopf, D. Vion, and F.K. Wilhelm for helpful discus- 
sions. This work was supported by the Air Force Office of 
Scientific Research under Grant F49-620-02- 1-0295, the 
Army Research Office under Grant DAAD-19-02-1-0187, 
the National Science Foundation under Grant EIA-020- 
5641, and the Advanced Research and Development Ac- 
tivity. 



[1] Y. Nakamura, Y. A. Pashkin, and J. S. Tsai, Nature 398, 

786 (1999). 
[2] D. Vion et al.. Science 296, 286 (2002). 
[3] J. Friedman et al.. Nature 46, 43 (2000). 
[4] C. H. van der Wal et al., Science 290, 773 (2000). 
[5] J. M. Martinis et al., Phys. Rev. Lett. 89, 117901 (2002). 
[6] R. C. Ramos et al., IEEE Trans. Appl. Supercon. 11, 998 

(2001). 

[7] J. E. Mooij et al.. Science 285, 1036 (1999). 
[8] B. L. T. Plourde et al., Phys. Rev. B 70, 140501(R) 
(2004). 

[9] I. Chiorescu et al.. Science 299, 1869 (2003). 
[10] I. Chiorescu et al.. Nature 431, 159 (2004). 
[11] Y. Makhlin, G. Schon, and A. Shnirman, Rev. Mod. 

Phys. 73, 357 (2001). 
[12] S. Saito et al., Phys. Rev. Lett. 93, 037001 (2004). 



5 



[13] E. H'ichcv et al, Phys. Rev. Lett. 91, 097906 (2003). 
[14] K. Lang et al., IEEE Trans. Appl. Supcrcond. 13, 989 
(2003). 

[15] J. M. Martinis, M. H. Devoret, and J. Clarke, Phys. Rev. 

B 35, 4682 (1987). 
[16] S. Linzen et al.. Rev. Sci. Inst. 75, 2541 (2004). 
[17] R. W. Simmonds et al., Phys. Rev. Lett. 93, 077003 



(2004). 

[18] P. A. Rcichardt ct al., Bull. Am. Phys. Soc. 49, A37.004 
(2004). 

[19] I. I. Rabi, Phys. Rev. 51, 652 (1937). 
[20] N. F. Ramsey, Phys. Rev. 78, 695 (1950). 



